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ABSTRACT 

We study the scattering problem of photon and polariton in a one-dimensional coupled-cavity system. Analytical approximate 
analysis and numerical simulation show that a photon can stimulate the photon emission from a polariton through polariton- 
photon collisions. This observation opens the possibility of photon-stimulated transition from insulating to radiative phase in 
a coupled-cavity QED system. Inversely, we also find that a polariton can be generated by a two-photon Raman scattering 
process. This paves the way towards single photon storage by the aid of atom-cavity interaction. 


Introduction 

A coupled-cavity QED system provides a promising platform to study novel quantum phenomena, since it combines two 
or more distinct quantum components, exhibiting features not seen in these individual systems. The discrete spatial mode 
of the photon in a coupled-cavity array and its nonlinear coupling to atom make the possible applications both in quantum 
information processing 1 and quantum simulation . 2 The seminal papers 3- - 1 proposed the use of the system to create strongly 
correlated many-body models. It has predicted the quantum phase transition from Mott insulator phase to superfluid phase . 5,6 
This scenario is constructed under the assumption that there is no extra photon leaking into the system. The stability of an 
insulating phase bases on the fact that the polariton states in a cavity QED system are eigenstates, i.e., spontaneous photon 
emission is forbidden. This situation may change if a photon can stimulate the photon emission from a polariton. In contrast 
with quantum phase transition induced by varying system parameters, such as atom-cavity coupling strength, stimulated 
photon emission from polaritons can also trigger the transition between insulating and radiative phases. It is interesting 
and important to investigate the photon-photon and photon-polariton scattering processes. Many efforts related to few-body 
dynamics mainly focused on multi-photon transports through coupled-cavity QED systems 7-19 , while a few works dealt with 
the formation of bound state . 20-22 So far, what happens when a photon collides with a polariton is still an open question. 

In this paper, we study the scattering problem of an incident photon by a polariton in a one-dimensional coupled-cavity 
QED system. Analytical approximate analysis and numerical simulation reveal several dynamical features. We find that 
a photon can stimulate the photon emission from a polariton, which induces the amplification of the photon population in 
a multi-polariton system. After a chain reaction, incident photons can stimulate the transition from insulating to radiative 
phases in the system with low doped cavity density. We also investigate the inverse process of stimulated photon emission 
from a polariton. We will show that a polariton can be generated by a two-photon Raman scattering process, which has been 
studied for the atoms found in nature . 23-25 Moreover, it has been shown that an atom-cavity system can behave as a quantum 
switch for the coherent transport of a single photon . 26 Considering a two-excitation problem, we find that a single-photon 
transmission through a quantum switch is affected significantly by a polariton that resides at it. 

This paper is organized as follows. At first, we present the model and single-excitation polaritonic states. Then, we 
propose an effective Hamiltonian to analyze the possibility of photon emission from two aspects. Numerical simulations for 
two-particle collision processes are showed later. Finally, we give a summary and discussion. 


Results 

Model and polariton 

We consider a one-dimensional coupled-cavity system with a two-level atom, which is embedded in the center of cavity array. 
The Hamiltonian can be written as 

N 

H = -KY j aja, +l +?ial\g)(e\ + H.c., ( 1 ) 
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where A represents atom-cavity coupling strength and K is the photon hopping strength for the tunneling between adjacent 
cavities. Here, |g) (|e)) denotes the ground (excited) state of the qubit with a z \e) = \e) and O z |g) = — |g), a/ (a} ) annihilates 
(creates) a photon at the /th cavity. Obviously the total excitation number, ,/K = L[/| _q a / a / + <J Z + is a conserved quantity 
for the Hamiltonian H, i.e., [H, - -V] = 0. 

The coupled-cavity array can be considered as a one-dimensional waveguide, while the two-level atom can act as a quan¬ 
tum switch to control the single-photon transmission. 26 To demonstrate this point, we rewrite the Hamiltonian in the form 

H = -2k£cos ka\a k + —== Y{ a l I#) (e| +H.c.) , (2) 
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where 
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It indicates that the atom couples to photons of all modes k G [—n,n\. In the -A = 1 subspace, atom can be regarded as a 
stationary scattering center. All the dynamics can be treated in the context of single-particle scattering method, which has 
been well studied. 26 

A comprehensive understanding for the dynamics involving the sector with -.■¥ > 1 is necessary to both theoretical ex¬ 
plorations and practical applications. Intuitively, the state of the atom (|e) or g)) should affect the interaction between the 
atom and a photon. In experiments, the practical processes may concern two or more photons, which obviously affect on the 
function of the quantum switch. On the other hand, the stability of an insulating phase may be spoiled by the background 
photons from environment. In this paper, we study the scattering problem in the jY = 2 sector, focusing on the effect of the 
nonlinearity arising from the atom. The investigation has two aspects: First, we study the photon scattering from a polariton. 
Secondly, we consider the collision of two photons under the atom-cavity nonlinear interaction. 

We start our investigation with the solution of single-particle bound and scattering states. In the invariant subspace with 
jY =1, exact solution shows that there are two bound states, termed as single-excitation polaritonic states, being the mixture 
of photonic and atomic excitations. From the Method, these polaritonic states are obtained by Bethe Ansatz method as the 
form 


^) = ±-^= sinh ^ ^ 1°^ + £ ^ r e ls> 1°) > 
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where the normalization factor is 
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a = 


2k\ z 
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The corresponding energy is 


£± = ±2K'coshj3, 

where the positive number ft determines the extension of bound states around the doped cavity, obeys the equation 
e 2p = \j(l/V2K^ +\ + (A / V2k) ‘ . 

We can see that j3 has nonzero solutions for nonzero A, indicating the existence of nontrivial bound states. 
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Figure 1 . (Color online) (a) Schematic configuration for the coherent collision of polariton and photon. An array of coupled 
single-mode cavities, where the central cavity is coupled to a two-level atom. Initially a polariton is located at the center, 
while a photon wave packet is moving from the left to collide with the polariton. (b) Schematic illustration for the equivalent 
description of the hybrid system. The excited state of the atom can be treated as a side-coupling site with infinite on-site 
repulsion. 


On the other hand, the derivation in Method shows that the solution of scattering states ]<//) with energy = —2ffcosk 
has the form 
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where A* is the normalization factor and 
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We can see that a polariton is a local eigen state of the system, which is stable and cannot emit a photon in the JV — \ 
subspace. The aim of this work is considering the effects of photon-photon and photon-polariton collisions. Our strategy is 
sketched in Fig. 1(a). In the invariant subspace with ,/K =2, a two-excitation state can be a direct product of a local photon 
and a polariton states, which are well separated in real space. As long as time evolution, two local particles are overlapped. 
The nonlinear effect induces the interaction between the photon and polariton. After a relaxation time, the free photons 
spread out from the central cavity, only the polaritons are left, being stationary at the center. In the case of the ultimate 
polaritonic probability being less than 1, (or the escaped photon number larger than 1) we can conclude that the stimulated 
photon emission occurs during the process. We will show that this behavior becomes crucial when we study the stability of a 
macroscopic insulating phase, and the efficiency of a quantum switch in a waveguide. In the following section, we will analyze 
the possibility of photon emission from two aspects. 


Effective description 

In this section, we present an analytical analysis on the effects of photon-photon and photon-polariton collisions. This will be 
based on an effective description of the original Hamiltonian H or H. We extend the Hilbert space by introducing the auxiliary 
photon state (aj)" |0) g , where a\ is the creation operator of a photon at site e and |0) e is the corresponding vacuum state. The 
qubit state \e) is replaced by a] |0) e . We rewrite the original Hamiltonians H and H as the Hubbard models 


H e q = — K E ajai + 1 + Xa^ae + H.c. + — a]a e (l — a\a e ), 
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and 


H e q = 



^E(«Ils)W + H.c.) 


—a\a e (l — a[a e ). 


(13) 


We note that the state (a|)" |0) e with n> 1 will be ruled out as U —> oo, Hamiltonians H eq and // cq being equivalent to H and 
H, respectively. Correspondingly, we have [H e q,od£ q ] = [77 e q,^eq] = 0 by defining ^ q =Lm =0 a / a / +a],a e . We will see that 
this equivalence can be true for a large magnitude U ~ 10. Next, we will perform our analysis from two aspects: k space and 
real space. 
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Coupled equations in k space 

First of all, we would like to point out that the eigenstates of the Hamiltionians H and H in Eqs. (1) and (2) are still the 
eigenstates of H eq and H eq by taking \e) —> a\ |0) e . Now we consider the case in two-particle subspace. The basis set for 
two-particle invariant subspace can be constructed by the single-particle eigen states 0 h ) and | </>*). We concern the complete 
basis set with even parity, which can be classified into four groups 


{|l,C7,k)} 
{|2,*,*')} 
{|3,ct,ct')} 



(14) 

(15) 

(16) 


where <7 = =b. We note that state | <j> k } (\) k j is automatically the eigenstate of 

| (f> a ) (j> a ^ will be ruled out as U —> ». Then basis sets (|l,a\k)} and {\2,k,k')} 
approximately. In this sense, the solution of the Schrodinger equation 


H with eigen energy £ k + £ k i. And states 
can further construct an invariant subspace 


ij- t \ ¥ (t))=H\y(t)), 


has the form 
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where coefficients C\ a k (t) and Ci/,-,/.-' (?) describe the two-particle dynamics and satisfy the coupled differential equations 
/—Ci(f) = M\\C\ (t) +M 12 C 2 (?), (19) 

«^C 2 (f) = M 22 C 2 (?) +M2\C\ (?). (20) 

Here the column vectors Ci (t) = {Cp^ (f)} and C 2 (?) = {Cj.k^' (?)}, and the matrix 


Mu M12 

M21 M22 


( 21 ) 


is a matrix representation of H on the basis set {Ci (t) ,C 2 (?)}■ Although we cannot get an analytical solution of |y/(f)}, 
we can conclude that the nontrivial solution | </?(?)) should predict the following relations in principle. We can always have 
nonzero C 2 (?) from initial condition Ci (0) ^ 0 but C 2 (0) = 0, i.e., 


in 


) 



and inversely, nonzero Ci (/) from initial condition C 2 (0) / 0 but C 1 (0) = 0, i.e.. 


( 22 ) 




(23) 


The former corresponds to the stimulated photon emission of the polariton, while the latter corresponds to the polariton state 
generation by a two-photon Raman scattering. The two processes are schematically illustrated in Figs. 2(a) and 3(a). 


Effective photon blockade 

In this section, we will demonstrate the process in Eq. (22) from an alternative way. One can consider the collision between 
an incident photon and an initial bound state around the site e in the system H eq . The obtained result should be close to that 
of the H system. In this context, the photon-photon collision only occurs at site e. Then the impact of the incident photon on 
the bound photon can be approximately regarded as a kicked potential on the <?th site. In the following, we will investigate the 
effect the potential works on the dynamics of the bound photon. 

We reduce the two-particle system of H eq to a single-particle system with the effective time-dependent Hamiltonian, 


4/14 



















Photon 


Polariton 


GS atom 


Photons 


© AAA/ 


(a) 


o AAA/ 
# AAA/ 


Insulating phase 


(b) 


?0 0 0 0 
0 0 o o o 


Radiative phase 

oO^o e>/or 

OjQ^O 


Figure 2. (Color online) Polariton-photon transition in a coupled-cavity array coupled to a two-level atomic system, (a) 
When the collision between a photon and polariton occurs, the total photon probability cannot be preserved. The gain of 
photons indicates the stimulated photon emission. The blue (gray) color represents the polaritonic (atomic ground) state, (b) 
The insulating-radiative phase transition. A multi-polaritonic insulating state can collapse to a radiative state by an external 
field radiation. 
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Figure 3. (Color online) The two-photon Raman transition in a coupled-cavity array coupled to a two-level atomic system, 
(a) When the collision between two photons from opposite directions occurs at the cavity with an atom, the total photon 
probability cannot be preserved. The loss of photons indicates the two-photon Raman transition. The blue (gray) color 
represents the polaritonic (atomic ground) state, (b) Single-photon storage by the aid of single photons train from the 
opposite side. Any single photons cannot be stored in the atom when they transmit unidirectionally. It can be achieved by an 
incident single photons from the opposite side. 
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H e ff (t ) 


(24) 


= H 0 + V(t), 

N 

H 0 = -k£ |/)</+l|+A|0)<e|+H.c., (25) 

|/|=o 

V(t) = U 0 8{t- z)\e){e\. (26) 

where Uq is the strength of the scattering and \e) = a\ |0) e , |/) = aj |0) ( l = 0, ±1, ±2,...) denotes the single-photon state. The 
initial state is one of the bound states 

|^) = -=[±- 5 - sinh/3 \e) + £ I/}]- (27) 

After the impact of the kicked potential, 10*} should probably jump to the scattering states <t) k ). In the following, we 
demonstrate this point based on time-dependent perturbation theory. 

For small Uq, the transition probability amplitude from the initial state |<p4) at t = 0 to |<p v ) (/i, v = ±) at t > T can express 
as 

A mv = — i J\<j) v \V (/') |</> M >e“*( e ^“ ev ) r, dr / (28) 

- E f dt ' f dt" e-^-^y ( 0 v |V (f') I 0 17 ) 

n ~k,± J o 70 

up to second order according to the time-dependent perturbation theory. Using the identity 

(^\V(t)\^) = (-I) 1 * 5 *- (0 J? |V(/)|0 v ), (29) 

and the completeness condition 


E <e|fWk) = l, (30) 

jj=&,± 

we get the transition probability between two bound states 

T+- = |A+_:| 2 = UqP 2 (l + Uq) , (31) 

T±± = |A±±| 2 = (l - Uqp) 2 + {Uop) 2 , (32) 

where 

4^2 

P = \e) (e\ 0 V > = (“I) sinh 2 (3- (33) 

The crucial conclusion is that the transition probability from the bound state to the scattering state is 

1 - r±± - r+_ = 21 / 0 / 7(1 - p - ufo) , (34) 


which is always positive for small nonzero Uq. This indicates that the collision between a photon and a polariton can induce 
the photon emission from the polariton. 

We employ the numerical simulation for verification and demonstration of our analysis. We compute the time evolution 
of an initial bound state by taking a rectangular approximation to a delta function. 


y(0 = 


^k)(e|, w > t T > 0 

0 , otherwise 


(35) 


For fixed Uq, we carry the calculation for different values of w. It is found that the result becomes stable as w decreases. 
The convergent data are adopted as an approximate numerical result. The evolution of an initial bound state under the central 


6/14 




I 


Figure 4. (Color online) Time evolution of the initial bound state |<t>(0)) = 0 ) in Eq. (27) driven by the 
rectangular-pulsed potential. The magnitude distributions of the evolved wave function \/P (J.t) for several instants are 
obtained as converging results for A = 0.8, Uq = 2, and w = 2 x 10 5 K' 1 . The red line indicates the probability of the <?th 
site. The circle(black) represents the initial profile of the bound state at site e and the circles(gray) represents the initial 
profile of the bound state as comparison to the profile of the final state. It shows that there is a particle probability spreading 
out from the center to the infinity, and the final bound state has almost the same shape as the initial one but less probability. 
This indicates that a kicked potential can induce the transition from the bound states to the scattering states. 


potential pulse is computed as well. The magnitude distribution of the evolved wave function \JP(l,t) = \ (I |<t>(f)}| is plotted 
in Fig. 4. Here the propose of using \JP {1,1) rather than the probability P (1.1) is to highlight the escaping wave packets from 
the center. We can see that there are two sub-wave packets propagating to the leftmost and rightmost, and the amplitude of 
the central bound state is reduced after this process. It can be predicted that the bound-state probability will keep decreasing 
by the successive pulses potential. 

The result of this section cannot be regarded as sufficient proof of the occurrence of the stimulated photon emission from 
a polariton. Nevertheless, it shows that there is a high possibility that such a process can happen. In the following section, we 
will investigate this phenomenon by numerical simulation. 

Numerical simulation 

In principle, one can explore the problem by solving the coupled equations (19) numerically. The truncation approximation 
is necessary since a numerous number of equations are involved. However, we can take an alternative way for truncation 
approximation, which is more efficiency for a discrete system. We can solve the Schrodinger Eq. (17) in finite real space by 
computing the time evolution of the initial state 


l^(0)) = |<P)|^ ), 


(36) 


where \<p) denotes local photonic state which is separated from polariton 0 ) in real space. The following analysis is also 
available for the state |0 + ). At time f, the evolved state is 


l®(0> = 


dIk' 

k,k! 


f 




k,n=± 




(37) 
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Figure 5. (Color online) Collision process between an incident photon wave packet and a polariton. The probability 
distributions ^ (Z,f) for several instants are obtained by the time evolution under the systems of (a) the Hamiltonian // eq with 
A = 2, U = 0 and (b) U = 10, (c) the original Hamiltonian H with A = 2 (or equivalently, // cq with U = °°). The incident 
wave packet has ko = n/2 and a = 0.3. The blue (red) dotted line indicates the initial probability of the 0th (<?th) site as 
comparison. It shows that the probability of the scattering photon is conserved for the non-interacting case with U = 0, but 
not conserved in the presence of nonlinearity in H. The result demonstrates the occurrence of stimulated photon emission 
from a polariton. Moreover, it is observed that the incident wave packet is totally reflected from the center in the case (a), but 
transmitted in the aid of the polariton. The very close similarity between (b) and (c) indicates that equivalence between 
Hubbard and the cavity-atom models. In both cases, the scattered and emitted photons are still local, keeping the similar 
shape of the incident one. 



Figure 6. (Color online) Emission probability from polariton sized Iq = 9, stimulated by the photon wave packet with 
a = 0.3 and different ko, for the system with A = 2. It shows that the transition probability can reach 0.4 at ko « 0.13k. 
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Figure 7. (Color online) (a) Schematic illustration for the scattering process of a moving wave packet and a stationary 
polariton at center of finite chain, (b) Plots of 1 — P ms ( t ) for the cases with Iq = 9, a = 0.3, A = 2, L = 120, ko = 3 tt/ 4 
(blue), it/2 (red). One can see that the probability converges to a nonzero constant at long-time scale. 


where \Q) denotes two-excitation polaritonic state. We consider the local photonic state \(p) as a Gaussian wave packet with 
momentum k o and initial center N A , which has the form 

I (P (N A ,k 0 )) = -^=Y. e ~^ r[, ~ NA?eika,a ) 1°) (38) 

v^o i 

where Oq = ®~F _i VA) 2 j s the normalization factor and the half-width of the wave packet is 2\f\n 2/ a. We take 2 VI n 2 /a <C 

1 

|A^ 4 1 to ensure the two particles being well separated initially. The evolved wave function |<t>(t)) is computed by exact 
numerical diagonalization. 

The probability distribution 

= (39) 

is plotted in Fig. 5 to show the profile of the evolved wave function. One can notice that in the photon-polariton collision 
process, the probability of the polariton is not conserved. This result has implications in two aspects: First, we achieve a 
better understanding of the occurrence of stimulated photon emission from a polariton. We find that the scattered and emitted 
photons are still local. This is crucial for the multi-polariton system, since the outcome photons can stimulate the photon 
emission of another polariton with high probability. Second, it provides evidence to support the equivalence between H eq with 
large U and the original H. 

The above result is for an incident wave packet with ko = n/2. We are interested in the dependence of emission proba¬ 
bility on the central momentum ko of the incident wave packet. The probability of the survival polaritons can be measured 
approximately by the photon probability within the region of the initial polariton resides in, i.e., 

h 

Fres (t) = £ (40) 

Ifeo 


where Iq denotes the extent of the polariton. Obviously, P res (t ) contains the probabilities of the residual polariton and the free 
photons within [— Iq, Iq], For infinite chain system, 1 — P KS (°°) equals to the photon emission probability T. In the numerical 
simulation, the system is finite, we take T = 1—Min^s (/)] within a finite time interval in order to avoid the error from the 
reflected photons. Results of F as function of ko presented in Fig. 6, show that the maximal photon emission probability 
reaches 0.4 at ko m 0.7371. We can see that the stimulated transition is significant, which indicates that a polariton is fragile 
against an incident photon. 

Now we explore a system with a portion of cavities with doped atom. For a well prepared insulating phase, which is 
formed by many independent polaritons, decreasing A can lead to the delocalization of the photons. The above analysis offers 
an alternative probability: external radiation can trigger a sudden change of the state. After the collision of an incident photon 
and the first polariton, the scattered and emitted photons can further stimulate other polaritons. In order to mimic such a chain 
reaction, we study the multi-collision process by computing the time evolution of the two-particle system in a long time scale. 
We consider a finite system, in which the scattered and emitted photons are reflected due to the open boundary condition. 
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It can simulate the repeating collision process, resulting in the continuous probability decay of polaritons. Results of our 
numerical simulations of 1 — P res (t) is presented in Fig. 7(b). It appears that the local average of P les (f) continuously decays 
at beginning as predicted and then converges to a nonzero constant. As pointed above, P les (f) may contain photon probability, 
leading to P les (t) > 1. However, the local maxima of 1 — P ms (t) can measure the stimulated transition approximately. 

We presume that a polariton should be washed out by successive collision. However, numerical result shows that the 
residual polariton probability does not tend to zero after a long time. There are two main reasons: First, as time goes on, any 
wave packets will spread, reducing the impact of photons on the polariton. Second, the inverse process of photon emission 
should be considered, in which two colliding photons can create polaritons. To demonstrate such a process, we compute the 
corresponding simulation. In this process, according to Eq. 36, the initial state can be expressed as 

1*^(0)} = |<p(^Va,?t/3)) \(p(—Na,—ti/3)) , (41) 


which implies that there are only two symmetry Gaussian wave packets at the beginning. At time t, the evolved state is 


l*(0> = e ,Ht 1^(0)) 


— 52 dkk' 

k,k? 


r 


4> k 


£,jU=± 




(42) 


where \t,) denotes the two-excitation polaritonic state. The probability distribution PS (l. t) at several typical instants is plotted 
in Fig. 8. One can see that in the photon-photon collision process, the probability of photons is not conserved as well, which 
indicates that a polariton can be created when two photons meet at the 0-th cavity. This shows that a polariton can be generated 
by two-photon Raman scattering. As a summary of numerical results, we conclude that a polariton cannot completely transmit 
to a photon by the collision from a single photon, and inversely, a photon cannot completely transmit to a polariton by the 
collision from a single photon. The essential reason is the energy conservation: two-photon energy cannot match that of one 
photon plus one polariton, i.e.. 


—2k cos k — 2k cos k 1 ^ — 2Kcosk" + e±. 


(43) 


This feature can also be employed to realize all-optical control of photon storage. One main task of quantum information 
science is to find physical implementations in which a flying qubit can be stopped to store or process quantum information. 
It has been shown that a flying qubit can be stopped and stored as a collective polariton by tuning the cavity-atom coupling 
strength adiabatically. 27 In the present cavity QED system, a single-photon wave packet can be a flying qubit, while a polariton 
can be regarded as a stopped photon, or a stationary qubit. Our result indicates that a single-photon wave packet, or a train 
of separated wave packets cannot excite a polariton if the atom is in ground state at the beginning. Then any incident single 
photons from one side cannot create a polariton solely, leaving the atom in the ground state. This can be expressed as equation 

( < l> ± \e- m fl\(p(Ni,k 0 ))=0, (44) 

i 

where A,- < 0, ko £ (0, n) and |A,+i — Ni\ 2\/ln2 /a, i.e., all the n wave packets incident from left and the neighboring wave 
packets are well separated. In contrast a photon can be stopped at the polariton with the aid of single photons train from the 
opposite side. This can be expressed as equation 

{^\e- m \(p (|A/b| ,-*o)>nip(tfb*o)> ^ °’ < 45 > 

i 

i.e, the atom partially absorbs a photon to form a polariton. The processes expressed by two above Eqs. are schematically 
illustrated in Fig. 3(b). 


Discussion 

In this paper, the scattering problem of photon and polariton in a one-dimensional coupled-cavity system has been theoretically 
investigated. The analysis shows that, a photon can stimulate the photon emission from a polariton, which suggests that the 
insulating phase is fragile against the external radiation for a system with a lower density of doped cavity. This result can 
have some applications in practice. For example, this provides a way to induce the amplification of the photon population 
in a multi-polariton system as a photon amplifier. On the other hand, we also find that two-photon Raman transition can 
occur in this cavity QED system, i.e., a stationary single-excitation polariton can be generated by three-body, two photons 
and atom, collision. This phenomenon can be used to design a scheme to stop and store a single photon. Although this two 
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Figure 8. (Color online) Collision process between two incident photon wave packets from leftmost and rightmost, 
respectively. The probability distributions £? (Z,f) for several instants are obtained by the time evolution under the system of 
the original Hamiltonian H with X = 2 (or equivalently, H eq with U = °°). The red line indicates the probability of the eth 
site. It shows that the probability of the scattering photon is not conserved in the presence of nonlinearity in H. It 
demonstrates the polariton can be created by the collision of two free photons. 


photon-polariton transitions is probabilistic, it reveals the peculiar features of two-excitation dynamics, which significantly 
differs from a single-particle scattering problem and opens a possibility to achieve all-optical control of a single photon. The 
underlying physics can be understood as the effective interaction of two photons arising from the nonlinearity in the doped 
cavity. These photon emission and absorption processes is an exclusive signature of correlated photons and could be applied 
to the quantum and optical device design. 


Methods 


The exact eigenstates with Jf = 1 

In this section, we present the exact eigenstates with ,jV = 1 for the Hamiltonian H. The Hamiltonian has parity symmetry 
[P,H] = 0, where Pa/P 1 = a_/. The odd-parity eigenstates can be obtained directly, which is 


r = 


~^= £ (sin kl)a] |g)|0) 

ijtQ 


V2N_f t t 

2 ,y TT k \ k a - 


Is) 10} 


with eigen energy fcy = — 2fceos£, where TT is the normalization factor and a ' k is the photon operator in k space, i.e., 
+ 1 £ 


|/|=o 


a] = 


V2N 

_Ly 

V2Nr 


e ikl aj, 


e~ ikl a f k . 


The solutions | (j) k ') with even parity are two folds: 
(i) For real k , the eigenstates has the form 


0 ") = gk \e) |0 )+f k al | g) |0) + £ (A^V +B k e-*Wa]) |g) |0), 


¥0 


where 


io)=n i°)/> fl /1°)/=°- 

l«l=o 


(46) 


(47) 

(48) 


(49) 


(50) 
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Submitting <f> k ) to the Schrodinger equation 


H 




we get the equations for coefficients g k , f k , A k , and B k , 
e k = -K(e ik + e- ik ), 

&k (^ k e ik + B k e lk ^j = — K ^A k e 2lk + B k e 2ik + fk) , 
e kfk = ~2 k (A k e lk + B k e ,k ) + fg k , 

£kgk = ffk■ 


The eigenstates \<j) k ) are two folds: 

(i) For real A, a straightforward derivation leads to 


A k 


B k 


fk 

£k 


-gke ,k 

4/ic 2 A sin k 

~gk 

4/ K A sin A 

gke lk 

4/ie 2 A sin A 

gk 

4/kA sin A 

£k gk 

A ’ 

—2fceosA. 


2k 1 e k + (A 2 - el) (e k + Ke ik ) 
2K£ k e~ lk — (A 2 — e 2 ) , 

2K 2 e k + (A 2 - el) (e k + Ke*') 
2K£ k e ,k - (A 2 —e|)j , 


Then we have 



£t = 


^{k>|o> + f4l*>|o> + 3 -^ 


± 


(A 2 - el) 2Ke k e :fik =f (A 2 - el) 


^Cte^VWlO)}, 

//o 


where A k is the normalization factor, and e k = £ k = —2)Ccos A. These are extended states. 

(ii) There are two eigenstates with complex A which can be seen as two bound states. The boundary condition 



= 0, for / —> ±<», 


and real e k require 


A k = 0 1 k = ij5 or 7t + ifi 


(51) 

(52) 

(53) 

(54) 

(55) 


(56) 

(57) 

(58) 

(59) 

(60) 
(61) 


(62) 

(63) 


with real /3 > 0. A straightforward derivation leads to 


Bk = fk , 

(64) 

A 2 = K 2 (e~ 2 P-e 2 P), 

(65) 

e± = ±2fccosh/3. 

(66) 

e 2fi = \j[l/V2 k) +\+(x/V2k) . 

(67) 


Then we have 
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e 


( 68 ) 


2k 


^> ± ) = ±-^= sinh P \ e > l°) + E 


l*l=o 


(Tl) ; 




\g)\0), 


where the normalization factor is 


a = 



2 

sinh 2 ft + coth j3. 


(69) 
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